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ABSTRACT 



The propagation of fire generated smoke into a shipboard space has been 
computationally modeled using a commercial code generated by Computational Fluid 
Dynamics Research Corporation (CFDRC). This study was based on space 01-163-2-L 
of an Arleigh Burke Class Flight IIA Destroyer. However, with changes, the model can 
be reconfigured to represent other shipboard spaces. Multiple smoke scenarios are 
applied to the space. For all scenarios, the inlet used is forward water tight door. Smoke 
enters the upper half of the door, while air enters through the bottom half. The 
temperature of the inlet fluids is altered to observe its effect on propagation. In the last 
scenario, the floor temperature is isothermally held at 1200 K to simulate a fire in the 
space below. The results of this scenario shows that extreme temperatures of adjacent 
spaces has minimal effect on propagation. The overall goal of this study is to show how 
computational methods can be used to model propagation of smoke in shipboard spaces. 
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I. INTRODUCTION 



A. BACKGROUND 

In 1987, USS STARK (FFG-31) was struck by two Iraqi Exocet missiles while on 
patrol in the Arabian Gulf. Even though the missiles failed to detonate, the crew had a 
variety of obstacles to overcome in order to save the ship. The largest problem facing 
them was the huge amount of unbumed fuel that was producing enormously high 
temperatures, and causing thermal damage to the ship’s mid-section. See Figure 1. 
Decks and electrical cableways, which melt at 2800°F (181 1°K), were no match for the 
3000° F(1922° K) produced by the unburned fuel. In addition to the tremendous amount 
of heat, the crew was forced to fight the blaze through an enormous amount of smoke. 
This hampered the crew’s ability to even get near the fight the fire and challenged the 
Navy’s fire fighting technology. 

The dense smoke severely limited the men’s ability to investigate spaces 
and to maintain sufficient breathing support. Early on, smoke completely 
engulfed the main-deck passageway, and many below-deck spaces had to 
be abandoned because of smoke before the impact of heat and flames was 
felt. Contributing to this difficulty was the toxicity of the smoke produced 
by the burning PVC-jacketed electrical cables and antisweat pipe 
insulation installed throughout the ship. [Ref 1] 

While the USS STARK had a state of the art fire protection system, it proved to 
be no match for the extreme temperatures of a “weapon induced” fire. This extremely 
high temperature fire is the greatest danger to any ship, even the Navy’s newest Arleigh 
Burke class destroyer. 



1 



THIS PAGE LEFT INTENTIONALLY BLANK 



2 



BRIDGE 










' 

CIC 


\ 


/ARDROOr 


\ 


RADAR EQUIR 
|_ ROOM H 




GALLEY 


:po berthing 


CREW BERTH] 




MK 13 
MISSILE 
MAG. 








CREW BERTHING 































■ FIRE BOUNDARY 

□ SMOKE/HEAT DAMAGE 

■ FIRE DAMAGE 



Figure 1. USS STARK (FFG-31). View of Starboard Side Damage 
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B. DESCRIPTION OF THE PROBLEM 



The engineering profession has realized that it is more cost effective to model and 
test a computer generated design than to build and test the same design in a laboratory. 
With the recent advances in desktop computers, engineers are now able to accurately 
model three-dimensional heat transfer and fluid flow problems. These advances could be 
used to design future ships by modeling shipboard spaces in order to determine the 
hazards of smoke resident in each design. If problems are discovered, improvements to 
the design could be made and then tests could be computed again to determine if the flaw 
still exists. This could save both time and money! 

This strategy could and should be used for fire protection systems resident in 
shipboard spaces. Fire protection engineers have not drastically changed the way in 
which ships are designed since the Perry Class was commissioned. With the speed and 
memory capability of the latest generation of desktop computers, problems can be 
solved, while they are still in a computationally generated form. 

Since the Surface Combatant of the 21 st Century (DD-21) is being designed for a 
maximum of 95 personnel, an enormous amount of planning will be required in order to 
ensure that the ship’s personnel will be able to safely defend the ship against fire and 
smoke. After all, fighting a fire is the most manpower intensive evolution that U.S. 
Navy ships undergo. The results of this study will give engineers another tool to 
determine realistic scenarios, which can be used to support or reject the proposed 
manning requirements. 
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C. PREVIOUS WORK 



Kay and Morris [Ref 2] examined the development of smoke modeling. They 
determined that fire safety aboard naval vessels was based on engineering judgement and 
previous fire incidents. However, with the speed of modem processors, computer 
modeling is possible, which provides engineers with a range of possibilities for solving 
the same problem. One can compare and contrast these solutions using either zone 
modeling or field modeling in the form of computational fluid dynamics (CFD). Using 
CFD there are two approaches to take: deterministic or probabilistic modeling. As 
expected, the probabilistic model quantifies the probability of fire spread. While these 
results may satisfactorily provide a quick answer in accordance with the laws of 
probability, shipboard fires spread based on the laws of nature, which may or may not 
coincide with the laws of probability. A better approach is to develop a physically based 
deterministic model, which produces more realistic informat ion(temperature, velocity, 
pressure etc). This method can better be assimilated to real world scenarios providing 
engineers with a broader approach, which is better suited to solve fire safety problems. 

In 1996, Tabarra, Kenrick and Matthews [Ref 3] analyzed natural smoke 
movement in a model tunnel using a CFD validation. They found “that a strong 
ventilation system may not be the best answer to passenger safety, as this tends to mix the 
hot smoke layers with the cool fresh air drawn in by the fire.” Their nemesis was the 
relatively low computational power of computers, but their coarse models accurately 
represented the flow. Unfortunately, radiation and turbulence were neglected, because 
their grids were already maximized to the computers’ limits. 
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Chow [Ref 4] used a CFD package, entitled PHOENICS, to computationally 
model a fire inside a building and how the resulting stratified smoke layer changed with 
different forced ventilation schemes. Also, using his models, he was able to predict the 
temperature of a well-mixed model if hot smoke is mixed with cool air. 

Null [Ref 5] used CFD-ACE to model the Concentric Canister Launcher(CCL) 
when it is exposed to a fire in an adjacent space aboard Navy ships. His research was 
valuable to predicting the temperatures inside the CCL, because once a cook-off 
temperature is reached, the energetic material (solid propellant, liquid fuels and 
explosives) in the missile could explode. The impetus behind his research was the 
unburned fuel of the Exocet missiles which hit the USS STARK(FFG-31) in 1987. The 
unburned fuel produces fires in excess of 3000° F (1922° K) 

To date, CFD research has focused on two-dimensional analysis, but with the 
recent advances in computational speeds and memory capacities, three-dimensional 
analysis is becoming feasible for researchers with just a desktop personal computer. 
Blay, Tuhault and Joubert [Ref 6] investigated a free convection problem in a corridor 
with smoke being removed through the ceiling via a soffit. The box in the lower left 
comer was a heat generating source, which was the driving force behind free convection. 
Their results included isothermal lines and velocity vectors of the smoke movement in 
the corridor. Although not stated, it is assumed that the CFD code, used to model the 
flow, was generated by the authors. 

D. OBJECTIVES 

The purpose of this study is to develop and examine a computationally generated 
model which will give engineers a better understanding of how smoke travels in 



7 



shipboard spaces. While the model will be generated in the form of passageway 01-163- 
2-L, onboard an Arleigh Burke (DDG-51) class destroyer, additional inlets and outlets for 
the smoke have been added to allow more opportunities to model other spaces or other 
scenarios. Furthermore, the source of the fire will determine the physical properties 
associated with the smoke, such as density, viscosity etc. These properties will have a 
dramatic effect on the movement of the smoke. 
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II. COMPUTATIONAL FLUID DYNAMICS 



A. OVERVIEW 

The software package used to analyze the smoke propagation was developed by 
CFD Research Corporation (CFDRC). It is entitled CFD-ACE [Refs 7-10] and uses three 
separate, yet interactive codes to solve the fluid flow problem. See Figure 2. The codes 
are GEOM, GUI and VIEW. CFD-GEOM is the pre-processor, which offers 
comprehensive geometry and grid generation, enabling the user to generate a grid 
representing the architecture of the problem. This stage will consume approximately 90 
percent of the total input time, due to the large number of inputs. CFD-GUI allows the 
user to input the specific thermophysical properties associated with the fluids. This 
includes the initial conditions (density, velocity), boundary conditions (adiabatic, 
isothermal) and interaction of species within the problem (heat transfer, turbulence, 
mixing). The user designates the number of iterations. Likewise, the results can be 
slightly varied by altering the differencing scheme, which will be explained in more 
detail later. Altering the relaxation method or amount of constraint can change the 
solution to the conservation equations. After all of the parameters have been entered, 
CFD-GUI solves the series of equations. The last stage, CFD- VIEW, is the 
postprocessor, which displays the results graphically through colored gradients, vectors 
or animation. 
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B. FINITE VOLUME ANALYSIS 



CFD-ACE divides the geometry into a finite number of volumes, each generated 
by a structured grid that uses single i, j, k coordinates to identify a particular direction 
within the finite volume. Various coordinate systems can be used in order to represent 
the problem ie. cartesian, cylindrical or spherical. Shipboard spaces are primarily 
rectangular, thus, a cartesian system was chosen. 

The following was taken from the Theory Manual of CFD-ACE and best 
describes the theory in the finite volume approach: 

This means that the discretized equations are formulated by evaluating and 
integrating fluxes across the faces of control volumes in order to satisfy the relevant 
conservation equations. Each control volume is, topologically, a cube with six faces and 
six direct neighbors. The cell centered aspect of the CFD-ACE approach means that each 
dependent variable is solved for at the center of each control volume and the values 
obtained are considered to prevail over the entire control volume. When evaluating 
convective fluxes over a control volume, differencing schemes of varying accuracy can 
be used. CFD-ACE offers a range of differencing schemes from first-order upwind, to 
second-order central differencing, to third-order schemes in high-speed flows. These 
schemes can be independently selected for each variable to be solved. 

Using a finite volume approach, the governing equations can be discretized into 
cells. CFD-ACE uses a cell centered variable arrangement where all flow variables and 
fluid properties are stored at the center of each cell. The following equation is derived 
for each variable in each finite volume: 

“ftp = a e4>e + “J* + a n<t>n + ° s<f> s + a h<t>h + a i<f>l + S 

The meaning of the symbols is as follows: ^ is the flow quantity, such as temperature, 
velocity or concentration. It is composed of average and instantaneous components. P 
refers to the center of the cell and e(East), w(West), n(North), s(South), h(high) and 
l(low) represent the values at the six adjacent faces of that cell, as shown in Figure 3. 
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The source term, S, represents terms other than the convection and diffusion effects, such 
as the pressure gradient. The coefficient ‘a’ represents the effects of convection and 
diffusion across the cell faces, which can be represented as follows: 



Convection Term - p e A e u e <j> e Diffusion Term 



-r 



0 ,- 6 ) 

A e 



K 



Ae is the area across the east face and A e is the distance from the center of the cell to the 
East face, p represents density, while Tis the molecular diffusion coefficient of the 
fluid. 



The dynamics of this fluid problem are driven by conduction, convection, mass 
diffusion and radiation. In order to understand the problem, the governing equations 
must be analyzed. Only the most general version of the equations will be used: 

r\n 

The Continuity Equation for compressible fluid is: — + div(pu) = 0 

dt 



The Continuity Equation for incompressible fluid is: 



du dv dw 
dx dy dz 



= 0 



In the Momentum Equation, p is the static pressure, Tij is the viscous stress tensor and f; is 
the body force. 



d.d dp 3r ;> 

dt ox, ox, ox. 



For Newtonian fluids, Tjjcan be related to the velocity gradients through 
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where p. is the fluid dynamic viscosity and 5jj is the Kronecker delta. Substitution of the 
above two equations results in the Navier-Stokes Equation: 



d 



Sp 5 i du i du j 2 du 






In the Energy Equation, shown below, Jjj is the total (concentration-driven and 
temperature driven) diffusive mass flux for species i. Also, hj represents the enthalpy for 
species i, S a stands for additional sources due to surface reaction, radiation and liquid 
spray and qj is the j-component of the heat flux. 

W + ±^M = -& + & + „±. + t *L-± ( j A)+S ' 

' dx j ^ J ' dxj dt 1 dxj ,J dxj dx, ,J ' a 



Fourier’s Law is employed to model the conduction heat flux where K is the thermal 
conductivity. 




CFD-ACE employs an iterative solution scheme in which the assembled 
equations for each dependent variable are solved sequentially and repeatedly with the 
goal of improving the solution with each iteration. This reduces the errors to acceptably 
small values. The nonlinear and coupled nature of the equations makes it necessary to 
restrain or under-relax the iteration-to iteration changes in the variable in order to obtain 
convergence of the solution procedure. 

The differencing scheme determines how the cell face values are calculated. The 
differencing scheme used in this model is the first-order upwind scheme. This means that 
the value of <f> at any given point is dependent upon the value of <f> at the adjacent point, 
just upstream of the flow. 
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For each iteration, the pressure (P) and velocity (U) terms will equal the current 
value plus some unknown correction according to the following equation: 

P=p* + p’ 

U=u* + u’ 

where P represents the current value of pressure, while p* represents the previous value 
of pressure as P is undergoing the iterative process and p’ is the unknown correction that 
exists between iterations. The symbols and process for computing the velocity remain 
the same as for the pressure. 

When the correction approaches zero then we can be reasonably confident that the 
value of pressure and velocity are accurate for the given flow field. The values of these 
corrections are determined by using the momentum equation to develop a functional 
relationship between the pressure and velocity fields. The pressures are then substituted 
for the velocities into a discretized continuity equation and the simultaneous matrix of 
equations is solved over the entire region. The solution procedure is summarized in the 
figure below. Figure 4 shows that for each iteration, within a transient time step, a 
system of equations must be solved for each dependent variable. A benefit to using the 
structured grid geometry is that it produces a banded matrix of coefficients for each 
dependent variable. The solution methods used for each system of equations are 
Backward Euler and a version of Forward Differencing Technique. 

In order to ensure that the solution for each iteration does not diverge, under- 
relaxation must be applied to the dependent and auxiliary variables. Under-relaxation 
constrains the amount that each variable can change from one iteration to the next. The 
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dependent variables (u,v,w,k) are modified using an Inertial Factor. A linear under- 
relaxation technique is applied to the auxiliary variable(p,T,// , p). 




Figure 4. CFD-ACE Flow Solver Solution Procedure [Ref 7]. 



15 



THIS PAGE INTENTIONALLY LEFT BLANK 



16 



III. MODEL 



A. GEOMETRY 

1. Model Selection 

The parameters for this model were taken from passageway 01-163-2-L of an 
Arleigh Burke Class Destroyer (DDG-21). The dimensions of the model identically 
match the actual passageway. Figure 5 shows a plan view of the compartment as shown 
in the ship’s drawings. This passageway (p-way) was chosen, because it is 
geometrically simple and has a variety of openings for smoke intrusion. The model 
allows a variety of smoke propagation scenarios to be studied. For example, even 
though the model may have an entrance in a specific location, unless that door is 
activated as an inlet or outlet, the program will treat it as a wall. Therefore, the non- 
activated entrances and exits will have no effect on the results. For example, this model 
could be used to represent the manner in which smoke enters a space, either horizontally 
or vertically. The model could also could be used to determine the rate or path which 
the smoke exits the space. This could be accomplished through a desmoking procedure 
using either forced or natural means, or it could be used to determine the affects of 
natural circulation in the skullery. 

2. Grid Distribution and Model Generation 

The simulations were carried out using a Micron Client Pro Desktop computer, 
with 128 megabytes of Ram and an internal 12 gigabyte hard drive. The software used 
was CFD-ACE version 5.0, which was last updated in February, 1999. 



17 




FR 
‘ 181 



|— -3* AFT FR 172 * 

| EXEC OFFICER | 




Figure 5. Plan View of Space 01-163-2-L 
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Generating a finite volume mesh is a sequential process. The first step is to 
generate a skeletal model as shown in Figure 6. This rough sketch shows the large 
rectangular outer structure of the space and the location of the hatches. 

The next step in generating a finite volume is to superimpose a grid upon each of 
the lines in Figure 6. Each major edge was allowed a total of 30 grid points. This meant 
that each of the smaller edges needed to be edged proportionally, according to its length. 
This formed an intricate puzzle, continuously summing the edges to ensure the totals of 
opposite edges matched. For most edges, this was not difficult due to the symmetrical 
nature of the planes perpendicular to the y and z axes. However, juggling the edge points 
of adjacent walls was difficult. 

Once the edging process was complete, faces were generated by designating four 
connected, edges. In most cases, rectangles resulted from connecting the four edges. If 
desired, any of the faces could be designated as an inlet or an outlet. This provides a 
great deal of the flexibility in the model. There are a total of 119 faces, as shown in 
Figure 7. This model allows the investigation of a variety of scenarios. Circles are used 
to represent scuttles (round, vertical entrances or exits to a shipboard space). However, 
circles need to be generated in a series of four arcs, which must follow the same 
restrictions as the four sides of a rectangle. These four arcs are again grouped in order to 
form a face. 

Once all of the faces have been generated, a block is formed. A block can best be 
defined as a group of six adjacent faces. For example, the rectangular volume, shown in 
Figure 8, has sides located at X=Y=Z=0 and an opposite face at their respective 
maximum values. 
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Figure 6. Line Drawing of Model using CFD-ACE 
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Figure 7. CFD Representation of P-WAY 01-163-2-L of DDG-51 Class (Extra 
hatches added to increase the number of possible scenarios) 
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Figure 8. The Black Lines designate the Inlets and Outlets of the compartment. 
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Once all six sides are designated, the code will recognize this new entity as a 
completely closed structure, called a domain. If desired, multiple domains could be 
linked to model a series of shipboard spaces. 

B. THERMOPHYSICAL MODEL 

After the physical dimensions have been constructed in CFD-GEOM, the domain is 

transferred into the numerical solution generation portion of CFD-ACE, entitled, CFD- 

GUI. The type of fluid flow problem is then determined: incompressible flow, heat 

transfer, radiation etc. Based on the type of problem entered, GUI requires inputs (initial 

conditions, boundary conditions) in order to solve the conservation equations. Many of 

the inputs are constant, such as the molecular weight. For, example, air is 29 g/mol and 

smoke is modeled as a gas with a molecular weight of 46 g/mol. However, others 

properties will change as the problem evolves. For example, the gaseous viscosity is 

computed by using Sutherland’s Law as shown below: 

AT 15 

u = 

B + T 

A and B are determined by the gases and T represents the temperature as it changes 
during the iterative process. The Joint Army, Navy, NASA, and Air Force (JANNAF) 
Tables [Ref 7] are used to estimate the specific heat, enthalpy and entropy for a gaseous 
species. When the Ideal Gas Law is selected, density is calculated as shown below: 

(P + Pref) M 

P~ 

RT 

Prefis the reference pressure, M is the fluid molecular weight and R is the universal gas 
constant. T and p are the locally calculated values of relative temperature and pressure 
respectively. 
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Turbulence is modeled using the standard k-s model. It is a two equation model 
that employs partial differential equations to govern the transport of turbulent kinetic 
energy, k, and its dissipation rate, 8. CFD-GUI utilizes the Launder and Spalding model 
(1974). Further details can be found in Ref 7. 

Spatial differencing allows the governing equations to be discretized, so the 
advantages of the finite volume method can be utilized. The finite volume method 
recognizes that the any property belonging to a finite volume will be influenced by all six 
adjacent cells. However, the upwind differencing scheme takes this idea one step further. 
It adds a weighting factor to the properties associated with each of the adjacent cells. The 
greatest weight goes to the cell, which is directly upwind of a given cell. Likewise, 
reduced weights are assigned to cells, which are downwind of a given cell, because their 
ability to change the given cell will be reduced due to the flowing stream. All scenarios 
used the upwind differencing scheme. 

C. REQUIRED INPUTS 

Due to the complex nature of CFD terminology, a brief description of the required 
inputs seems appropriate. 

1. Under-relaxation 

Under-relaxation is used to constrain the change in the variable from iteration to 
iteration in order to prevent divergence of the solution procedure. Large values of under- 
relaxation are associated with variables that have an increased risk of becoming 
divergent. 
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2. Sweeps 

Sweeps can best be described as a local iteration within the global iteration 
process. As the global finite volume process iterates, the solver will iterate a given cell 
an increased number of times in order to prevent the solution from diverging. For 
example, the default setting for all velocity sweeps is 5. However, pressure can cause a 
solution to diverge quicker than other variables, therefore it’s default setting is 30. 

3. Initial and Boundary Conditions 

The grid is based on a cartesian coordinate system with U, V, and W velocities 
aligning themselves with the X, Y, and Z axes, respectively. All pressures in CFD-ACE 
are relative to a reference pressure, which is set at 100000 N/m 2 . Qualitatively, 
Turbulence Kinetic energy (K) is total amount of energy associated with an eddy. Just as 
the title indicates, Turbulence Dissipation Rate (D) is the rate at which the eddy loses 
energy. Turbulent Length Scale (L) is estimated at approximately 0.3% of the size of the 
inlet. Only two of the previous three variables are required to process the solution, while 
the third is computed using the equations below: 

11.74X 15 

L 

If all three variables are entered, L is ignored and generated via the above equation, 
using D and K. 
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IV. RESULTS 



The main objective of this research was to use a commercially-developed CFD 
program to generate a three-dimensional model of the smoke dissipation within the 
configuration of a shipboard space. 

With the assumption that this is the first research to be done coupling this CFD 
software with the goal of shipboard space modeling, the accuracy of the results with 
respect to a live testing is unclear. However, one internal method of measuring the 
accuracy is to view the residual plots. For each scenario completed, the residuals all 
decreased at least four orders of magnitude, which meets CFDRC’s criteria for 
convergence. Therefore, the results can be considered valid. 

Based on the success in the residual converging for each scenario, it leads this 
engineer to believe that the CFD code has not been pushed to its limits. Most scenarios 
in fluid flow require alterations of the solution monitors (sweeps or under-relaxation) in 
order for a solution to converge. This study regularly used the recommended values for 
sweeps and under-relaxation and the solution still converged. This indicates that a more 
violent flow scenario (higher velocities or increased turbulence) could have been 
modeled before the solution failed to converge. Of course, this would then require the 
inherent monitors to be altered to regain a converging solution. 

As with any modeling software, the results are only as good as the input data. 
CFD-ACE is no exception. The user needs to be aware that the slightest change in any 
input data can drastically change the output. Other dramatic changes can be viewed by 
changing the iterative solving mechanisms, such as the relaxation or sweep parameters. 
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When contrasted, scenarios A, B and C exhibit small differences in the diffusion 
rates. This is attributable due to the greater molecular mass of the smoke, 46 g/mol, 
compared to the air, at only 30 g/mol. 

Contrasting Scenarios A and C, it is apparent that the dissipation rate is not 
affected, even when the floor is exposed to such extreme temperatures. 
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V. CONCLUSIONS 



This study’s results successfully modeled smoke propagation in shipboard spaces. 
The accuracy of the results is valid as verified through the use of residuals. 

This study successfully pioneered the use of today’s CFD technology to simulate 
the effects of smoke dissipation in shipboard spaces. These models were generated in 
order to develop a greater appreciation for the numerous forces that alter smoke 
movement. As the power of computer and increased accuracy of CFD technology grow 
together, this method of shipboard simulation will one day become the norm for system 
design and engineering. 
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VI. RECOMMENDATIONS 



The following recommendations are made in the continuation of this study: 

Run the model using a fixed pressure at the inlet and the outlet. The Arleigh Burke 
Class has a unique CBR pressurizing system, which will most likely be the driving 
force behind smoke propagation. 

- Change the geometry to reflect the presence of the ladder, which runs between the 
upper and lower hatches. 

- Use CFD to generate fire scenarios, which will aid in determining the manning 
requirements for DD-21, which is still in the design phase 

Develop a scenario which has a fire within the modeled space to determine its effects. 
Add on adjacent spaces to see how the fire will affect them. 

Use CFD to compute the time it takes the wall to reach a given temperature using the 
transient function that will be available in version 6.0 
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APPENDIX A. 



Scenario A was developed as a control Scenario. The results of this scenario will 
be compared to Scenarios B and C. The inlet is located at the X-max face with the upper 
half of the door designated as the smoke inlet, while the lower half is exclusively air. 
Refer to the next page for the required inputs. 
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Relaxation 


Velocity 


(m/s) 


0.2 




Turbulence 


(J) 


1.0 




Enthalpy 


(KJ/kg) 


1.0 




Mixture Fractions 


0.2 


Sweeps 


Velocity 


(m/s) 


5 




Pressure 


(Pa) 


30 


Initial 

Conditions 


U Velocity 


(m/s) 


-0.1 




V Velocity 


(m/s) 


0 




W Velocity 


(m/s) 


0 




Relative pressure 


(Pa) 


0 




Turbulence Kinetic Energy 


(J) 


0.04 




Rate of Turbulence Dissipation 


(J/s) 


-0.05 




Turbulent Length Scale 


(m) 


0.06 




Temperature 


(K) 


500 




Reference Pressure 


(Pa) 


1E5 


Boundary 

Conditions 


Isothermal Wall Temperature 


(K) 


300 


Inlet - Smoke 


U Velocity 


(m/s) 


-0.1 




V Velocity 


(m/s) 


0 




W Velocity 


(m/s) 


0 




Temperature 


(K) 


500 




Turbulence Kinetic Energy 


(J) 


0 . 




Rate of Turbulence Dissipation 


(J/s) 


0.04 




Turbulence Length Scale 


(m) 


0.06 




Pressure 


(Pa) 


0 


Inlet - Air 


U Velocity 


(m/s) 


-0.1 




V Velocity 


(m/s) 


0 




W Velocity 


(m/s) 


0 




Temperature 


(K) 


500 




Turbulence Kinetic Energy 


(J) 


0.04 




Rate of Turbulence Dissipation 


(J/s) 


-0.05 




Turbulence Length Scale 


(m) 


0.06 




Pressure 


(Pa) 


0 



Table 1. Input Data for Scenario A ( See Chapter III for definitions) 
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Outlet 


U Velocity 


(m/s) 


0 




V Velocity 


(m/s) 


0 




W Velocity 


(m/s) 


0 




Temperature 


(K) 


400 




Turbulence Kinetic Energy 


(J) 


0.02 




Rate of Turbulence Dissipation 


(J/s) 


-0.05 




Turbulence Length Scale 


(m) 


0.08 




Pressure 


(Pa) 


0 



Table 1 Cont. Input Data for Scenario A ( See Chapter m for definitions) 




Figure 9 Three dimensional design of Inlet (Red door on Right) and Outlet (Red 
Door on Left). Top Half of Inlet is Smoke while Bottom Half is Air. 
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Figure 10. Isosurface of Smoke showing 98 percent concentration by volume 
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Figure 11. Isosurface of Smoke showing 54 percent concentration by volume 
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Figur e ,2. — ce of Smoke showing 22 percent b y vo.ume 
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Figure 13. Isosurface of Air showing 77 percent by volume 
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Figure 14. Isosurface of Air at 54 percent concentration by volume 
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Figure 15. Isosurface of Air at 16 percent concentration by volume 
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Figure 16. Residuals for Scenario A 
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APPENDIX B 



Scenario B was developed to contrast with Scenario A. The inlet is located at the 
X-max face with the upper half of the door designated as the smoke inlet, while the lower 
half is exclusively air. The major difference between this model and the others is that the 
smoke is 500 K, while the air is 300 K. Refer to the next page for the required inputs. 
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Relaxation 


Velocity 


(m/s) 


0.2 




Turbulence 


(J) 


1.0 




Enthalpy 


(KJ/kg) 


1.0 




Mixture Fractions 


0.2 


Sweeps 


Velocity 


(m/s) 


5 




Pressure 


(Pa) 


30 


Initial 

Conditions 


U Velocity 


(m/s) 


-0.1 




V Velocity 


(m/s) 


0 




W Velocity 


(m/s) 


0 




Relative pressure 


(Pa) 


0 




Turbulence Kinetic Energy 


(J) 


0.04 




Rate of Turbulence Dissipation 


(J/s) 


-0.05 




Turbulent Length Scale 


(m) 


0.06 




Temperature 


(K) 


400 




Reference Pressure 


(Pa) 


1E5 


Boundary 

Conditions 


Isothermal Wall Temperature 


(K) 


300 


Inlet - Smoke 


U Velocity 


(m/s) 


-0.1 




V Velocity 


(m/s) 


0 




W Velocity 


(m/s) 


0 




Temperature 


(K) 


500 




Turbulence Kinetic Energy 


(J) 


0.04 




Rate of Turbulence Dissipation 


(J/s) 


-0.05 




Turbulence Length Scale 


(m) 


0.06 




Pressure 


(Pa) 


0 


Inlet - Air 


U Velocity 


(m/s) 


-0.1 




V Velocity 


(m/s) 


0 




W Velocity 


(m/s) 


0 




Temperature 


(K) 


300 




Turbulence Kinetic Energy 


(J) 


0.04 




Rate of Turbulence Dissipation 


(J/s) 


-0.05 




Turbulence Length Scale 


(m) 


0.06 




Pressure 


(Pa) 


0 



Table 2. Input Data for Scenario B ( See Chapter III for definitions) 
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Outlet 


U Velocity 


(m/s) 


0 




V Velocity 


(m/s) 


0 




W Velocity 


(m/s) 


0 




Temperature 


(K) 


400 




Turbulence Kinetic Energy 


(J) 


0.02 




Rate of Turbulence Dissipation 


(J/s) 


-0.05 




Turbulence Length Scale 


(m) 


0.08 




Pressure 


(Pa) 


0 



Table 2 Cont. Input Data for Scenario B ( See Chapter III for definitions) 
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Figure 17. The Inlet is on the right with smoke (blue) as the upper fluid at 500 K 
and air (red) on the bottom at 300 K. The Outlet is shown on the left. 
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Figure 



18. Isosurface of Smoke showing 



98 percent concentration by volume 



53 



CH6N2 

1 




Figure 19. Isosurface of Smoke showing 54 percent concentration by volume 
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Figure 20. Isosurface of Smoke showing 22 percent concentration by volume 
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Figure 21. Isosurface of Air at 77 percent concentration by volume 
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Figure 22. Isosurface of Air at 54 percent concentration by volume 



57 



02 




Figure 23. Isosurface of Air at 16 percent concentration by volume 
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Figure 24. Residuals for Scenario B 
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APPENDIX C 



Scenario C was developed to contrast with Scenarios A and B. The inlet 
is located at the X-max face with the upper half of the door designated as the smoke inlet, 
while the lower half is exclusively air. Both fluids enter at 500 K. The Y-min face is 
held isothermally at 1200 K. Refer to the next page for the required inputs. 
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Relaxation 


Velocity 


(m/s) 


0.2 




Turbulence 


(J) 


1.0 




Enthalpy 


(KJ/kg) 


1.0 




Mixture Fractions 


0.2 


Sweeps 


Velocity 


(m/s) 


5 




Pressure 


(Pa) 


30 


Initial 

Conditions 


U Velocity 


(m/s) 


-0.1 




V Velocity 


(m/s) 


0 




W Velocity 


(m/s) 


0 




Relative pressure 


(Pa) 


0 




Turbulence Kinetic Energy 


(J) 


0.04 




Rate of Turbulence Dissipation 


(J/s) 


-0.05 




Turbulent Length Scale 


(m) 


0.06 




Temperature 


(K) 


500 




Reference Pressure 


(Pa) 


1E5 


Boundary 

Conditions 


Isothermal Floor Temperature 


(K) 


1200 




Isothermal Wall Temperature 


(K) 


300 


Inlet - Smoke 


U Velocity 


(m/s) 


-0.1 




V Velocity 


(m/s) 


0 




W Velocity 


(m/s) 


0 




Temperature 


(K) 


500 




Turbulence Kinetic Energy 


(J) 


0.04 




Rate of Turbulence Dissipation 


(J/s) 


-0.05 




Turbulence Length Scale 


(m) 


0.06 




Pressure 


(Pa) 


0 


Inlet - Air 


U Velocity 


(m/s) 


-0.1 




V Velocity 


(m/s) 


0 




W Velocity 


(m/s) 


0 




Temperature 


(K) 


500 




Turbulence Kinetic Energy 


(J) 


0.04 




Rate of Turbulence Dissipation 


(J/s) 


-0.05 




Turbulence Length Scale 


(m) 


0.06 




Pressure 


(Pa) 


0 



Table 3. Input Data for Scenario C ( See Chapter III for definitions) 
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Outlet 


U Velocity 


(m/s) 


0 




V Velocity 


(m/s) 


0 




W Velocity 


(m/s) 


0 




Temperature 


(K) 


400 




Turbulence Kinetic Energy 


(J) 


0.02 




Rate of Turbulence Dissipation 


(J/s) 


-0.05 




Turbulence Length Scale 


(m) 


0.08 




Pressure 


(Pa) 


0 



Table 3 Cont. Input Data for Scenario C ( See Chapter III for definitions) 
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Figure 25. The Inlet is on the right with smoke (blue) as the upper fluid at 500 K 
and air (red) on the bottom at 500 K. The Outlet is shown on the left. 
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Figure 26. Isosurface of Smoke showing 98 percent concentration by volume. 
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Figure 27. Isosurface of Smoke showing 54 percent concentration by volume. 
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Figure 28. Isosurface of Smoke showing 22 percent concentration by volume. 
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Figure 29. Isosurface of Air showing 77 percent concentration by volume. 
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Figure 30. Isosurface of Smoke showing 54 percent concentration by volume. 
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Figure 31. 



Isosurface of Smoke showing 16 percent concentration by volume. 
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Figure 32. Residuals for Scenario C 
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